Wigner function for the orientation state 
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We introduce a quantum phase space representation for the orientation state of extended quantum 
objects, using the Euler angles and their conjugate momenta as phase space coordinates. It exhibits 
the same properties as the standard Wigner function and thus provides an intuitive framework for 
discussing quantum effects and semiclassical approximations in the rotational motion. Examples 
illustrating the viability of this quasi-probability distribution include the phase space description of 
a molecular alignment effect. 
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Introduction. — The successful treatment of a physical 
problem often rests on the right choice of formalism, and 
conceptual progress is hard if one chooses an inappro- 
priate language. This is particularly evident in quantum 
theory, where at first glance diverging, yet equivalent for- 
malisms come along with inherent advantages and disad- 
vantages. For good reasons, the hydrogen atom is usually 
quantized with Schrodinger's wave equation, while fields 
are naturally described by Heisenberg operators. Feyn- 
man's path integral formulation, on the other hand, is 
often preferred in semiclassical considerations, and vast 
perturbative expansions are best handled diagrammati- 
cally. It is clearly helpful to command a variety of alter- 
native methods for treating a given physical system. 

Wigner's prescription pQ to transform the density oper- 
ator into a distribution on classical phase space initiated 
yet another equivalent description of quantum mechan- 
ics. In particular, no other formalism is as powerful in 
highlighting quantum features while providing us with an 
intuition in classical terms. Moyal completed it [2] by dis- 
covering that the one-to-one mapping between quantum 
operators and phase space functions reflects Weyl's cor- 
respondence rule [3]. Since then, this picture has found 
numerous applications in many areas of physics, ranging 
from solid state physics to quantum optics 0HS], from 
particle interferometry to molecular physics [T0HI2] . and 
from semiclassics to the foundations of quantum mechan- 
ics [I3HIE]. 

Given the success of the standard Wigner- Weyl repre- 
sentation for point particles, it is striking that a viable 
phase space description for spatially extended quantum 
objects is still unknown. Such an extension to the orien- 
tation state of a rigid body should be particularly useful 
for the growing number of experiments probing quantum 
and classical dynamics in the rotation and alignment of 
complex molecules. There are prescriptions to construct 
Wigner functions from underlying Lie groups fH)H2S] , but 
applied to the orientation state one is not able to recover 
all essential features of the standard Wigner function, be 
it the quasi-probability interpretation or the direct cor- 
respondence to the relevant classical phase space. 

In this letter we present a Wigner function for the ori- 



entation state of an extended quantum object. In con- 
trast to the approaches mentioned above, it is founded 
on the quantized canonical coordinates of the underlying 
classical phase space, providing a one-to-one mapping be- 
tween quantum operators and classical phase space func- 
tions, the Weyl symbols. In close analogy to the stan- 
dard Moyal formulation, we find that any Weyl-ordered 
observable is mapped to its equivalent function on phase 
space, that integrating out phase space coordinates yields 
the reduced probability distribution, and that the motion 
of the Wigner function is described by a quantum Liou- 
ville equation, which turns into its classical equivalent as 
0. 

Wigner function for the orientation. — The Hilbert 
space % of a rigid rotor is spanned by the vectors 
\JKM), the eigenstates of the square of the angular 
momentum operator J 2 , and its projections P z and J z 
on the body-fixed and space-fixed z-axis. An alterna- 
tive basis is provided by the resolution of the identity 
I = J 27r da Jq d/3sin/3 d7|cv, j3, 7) (a, j3, 7I in terms of 
the Euler angles (a, (3, 7) = f2, which serve to specify the 
orientation (in the usual z-y-z convention). The basis 

transformation (f2| J KM) = \J "^rD J ^ K (Vl) is mediated 

by the Wigner D-matrices D J MK (fl) [26] . 

Choosing the Euler angles as natural configuration 
space coordinates, the classical phase space is completed 
by their conjugate momenta p a ,pp,p 1 . As the genera- 
tors of the rotations defining the Euler angles they are 
projections of the angular momentum vector J onto the 
corresponding rotation axes: p a is the projection of J 
on the space-fixed z-axis, pp on the nodal line formed by 
intersecting the space- fixed x-y plane with the body-fixed 
x-y plane, and p 1 on the body-fixed z-axis. 

Quantum mechanically, the symmetrized projections 
are given by the operators p Q = —ihd a , pp — 
— ih [dp + I cot (3) , and p 7 — —ihdj 26 . Like J 2 , J z , 
P z , they form a complete set of commuting observables 
exhibiting a discrete spectrum, with eigenvalues p a — 
hm a ,pp — 2hmp and p 7 — hm 1 , m a ^ n € Z. We denote 
the corresponding eigenstates as |m ct m ) gm 7 ) = |m); they 
provide a further basis of % and their Euler angle repre- 
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sentation is (Q\m) = e im " a e am ^e im -it /y/^ sin/3. 

With this, we can now state our Wigner function 
W (fl,m) = W (a, (3, 'y,p a /h,pp/2h,p 1 /H) in terms of 



Here we use the abbreviations a'± = (a ± a' /2) mod 2ir, 
alike for 7, (3' ± = ((3 ± (3'/2) mod7r, and sine (a?) = 
njSin(xi) /xi. While the angular representation Q is 
similar to the standard Wigner function [T] , the momen- 
tum representation ^ differs notably exhibiting a dou- 
ble sum and sine functions, which act like blurred Kro- 
necker deltas. Note that the angular momentum argu- 
ments m are integers reflecting the discrete spectra of 
the corresponding operators. When discussing the semi- 
classical limit, and for graphical representations, it is ad- 
vantageous to extend the definition to real arguments 
by rounding them to the nearest integer. The Wigner 
function in terms of the symmetric top eigenstate basis 
(JKM\p\J'K'M') is readily obtained from |l]) by insert- 
ing the corresponding resolution of the identity. 

As in the standard case, one can introduce a 
Stratonovich-Weyl operator kernel A (r2,m) (see below) 
|27j . such that the mapping of any operator A onto its 
Weyl symbol Wa {Q, m) reads as 

W A (0, m) = tr [AA (Q, m)] (3) 

and its inversion 

A=— V/ da df3 d 1 W A (n,m)A(n,m). 
47r m J o Jo Jo 

For A = p this implies (1) and (2), since W p (Q.,m) = 
4ir 3 W(n,m). 

Relevant properties of the Wigner function. — 
Remarkably, all relevant properties of the stan- 
dard Wigner function can be recovered. It follows 
from ([I]) that the Weyl symbol of any hermitian 
operator is real, as required for a proper phase 
space representation. Moreover, W (il, m) is nor- 
malized with respect to the phase space integral, 
j 'dOX; m W (0,m) = 1, where dQ = dadfidj, and 
expectation values can be calculated by the phase space 
average fdClY, m W A (fi,m) W (Q,m) = (A). The 
quasi-probability interpretation of the Wigner function 
manifests itself in the expressions for the probability dis- 
tributions (f2|p|fi) and (m\p\m). They can be obtained 
as marginal distributions of W (fl, m) by integrating out 
the conjugate variable, Y^ m W (f2,m) = sin /3{Sl|yo|0) 



both the orientation representation (U!\p\Vt") and the an- 
gular momentum representation (m'\p\m") of the quan- 
tum state, 

(1) 
(2) 

I 

and J W (fl, m) dVl — (m\p\m), as is readily verified 
from ^ and 

As for the Weyl symbols of the coordinate and momen- 
tum operators, Eq. ([I]) yields the expected expressions 
W a (f2,m) = a and W Pa (f2,m) = hm a = p a . More- 
over, in close analogy to the standard case, the operator 
ordering {& n ,p™} w = 2"™ £™ =0 (™)p™~ k a n p h a turns 
the Weyl symbol of arbitrary products of a and p a into 
their classical equivalents, ^{q™^™}^ (fi,m) — a ™P™; 
the other angles satisfy analogous expressions. Unlike in 
the standard case, this 'Weyl' ordering is not equivalent 
to the 'symmetric' ordering due to the operator-valued 
commutator relation [a, p a ]. 

Construction of the Wigner function. — To 
shed light on how the Wigner function is de- 
vised, let us consider the displacement operators 
D a (a,m a ) — exp(im a a)exp(—iap a /h), similar for 
7, and Dp(j3,mp) — exp(i2m^/3) exp (— i/3pp/h), 
where angular values outside the Euler range 
are mapped back by taking the modulus, 
e.g. D a (a, 0) |a ) = | (a + a) mod27r). These phase 
space translation operators commute, as follows from 
the commutators of the phase space coordinates. The 
D (f2, m) = D a (a,m a ) ® Dp(/3,mp) <g> D 7 (7,to 7 ) 
allow one to construct the major building block of 
the phase space formulation, the operator kernel 
A(0,m) = D(n,m)A(0,0)Dt(O,m), where A (0,0) 
is the direct product of 

A « (°> °) = ^ E f da ' D « m ) e~ W/2 , (4) 

A /3 (0,0) = i^ I '' d/fDp (13', m)e- im P', (5) 

and A 7 (0, 0) having the same form as Q. The choice of 
the phase factor in Q, ^ guarantees that the angular- 
momentum symmetry of a state is correctly reflected by 
the Wigner function 28J . The symmetric choice of the in- 
tegral limits, on the other hand, is required to ensure the 
hermiticity of the A, (0,0). Note that the negative lower 
integration boundaries do not contradict the definition 
range of the Euler angles since they refer to translations. 



W(Sl,m) 



da' 



47T 3 

1 y 



7T/2 
-7T/2 



47T 3 



dp' I d 7 ' ^/sin p\ sin (3'_ e im » a ' e l2m ? ' e m ^' (a'_ , P'_ , i_ \ p\ a' + , f3' + , 7^) 
e»(<-<') Q e i2 ( m ! 3 ~ m ! 3 ')/ 3 e l ( m: '" m ")T'(m'|p|m''). 



m + m 
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Inserting the explicit form of A (fi, m) in ^ yields the 
expressions given in (0) and ^. 

These results reproduce a further crucial property of 
the Wigner function: it properly reflects phase space 
translations in the sense that the Weyl symbol of a trans- 
lated operator A' = D (ft' , m') AD^ (f2', m') is the trans- 
lated Weyl symbol Wa — Q', m — m'). This is proven 
by rewriting tr [A'A (fl, m)] as tr [AA (fi - Q',m - m')] 
using the invariance of the trace under cyclic permuta- 
tions. 

Quantum Liouville equation of a free symmetric top. — 
Applying the definition ([lj to the von Neumann equa- 
tion allows one to determine the quantum Liouville equa- 
tion, the law of motion for the Wigner function. To 
this end, the free symmetric top Hamiltonian [55] must 
be recast in terms of the canonical operators, H = 
(p Q - p 7 cos/3) 2 /(2/isin 2 $) + p2/2h + P 2 /2I 3 -h 2 {l + 
sin -2 $)/8Ii, with 1\ and 1$ the principal moments of 



inertia. One thus recovers the classical Hamiltonian with 
an additional quantum potential given by the last term. 
Owing to this quantum correction one does not recover 
the classical Liouville equation for the torque-free mo- 
tion, unlike the case of a free point particle. This is 
a consequence of quantization and not inherent to the 
chosen phase space approach. Moreover, the discrete- 
ness of the conjugate momenta yields summations over 
m, in variance with the form of the classical Liouville 
equation. However, in the limit that the Wigner func- 
tion varies only weakly in the discrete momenta they can 
be replaced by continuous variables. This is the case 
for a "classical" state with macroscopic extension over 
phase space. If one further rescales the dimensionless 
Wigner function to W (Q,p) = W (fi, m) / (2h 3 ) (where 
Pa,Pp,P-y replace hm a , 2hmp and Km-y), the quantum Li- 
ouville equation assumes the form of the classical Liou- 
ville equation to leading order in H, 



d t w(n, P ) = - 



h sm (3 da h op 

j— q— - (p~ - p a COS 0) (p a - p 1 COS (3) 

1\ sin p 



h sin 2 /3 

) — 
J dpp 



(p 7 COS 2 /3 — p a COS (3) + — 

w{n, P ) + o(h 2 ). 



d 

d-j 



(6) 



Hence, the classical dynamics is retained up to correc- 
tions of order ti 2 , like in the standard case of a point 
particle. 

Representative Wigner functions. — As a first 
application, we consider the Wigner function 
of a coherent state in a [30] . \a,m a ) a = 
D (a, m a ) |0, 0) a <& \mp) ® \m 7 ), where |0,0) Q = 

■d^ 1/2 (0,exp [-l/o- 2 ]) E ma ex P {- m a/ (J ' 2 ) \m a ) defines 
the coherent state at the phase space origin (with $3 
a Jacobi theta function). Figure [l] shows the coherent 
state |vr, 10) Q with momentum spread a = 7. As one 
expects, the Wigner function is well localized in a and 
m a . If one decreases the angular width the discreteness 
in nia gets even less pronounced, illustrating that m a 
can be replaced by a continuous variable in the classical 
limit. We note that, unlike the standard Wigner 
function for Glauber coherent states, the phase space 
function of |tt, 10) q takes on negative values in parts 
of the phase space. However, since these are negligibly 
small compared to the peak height in Fig. [l] one can 
maintain that the coherent states provide a classical 
correspondence. Coherent states involving more than a 
single angle exhibit a similar behavior. 

If one superposes two coherent states separated in 
phase space, oscillatory interference terms emerge in be- 
tween the classical contributions, see Fig. [2} We find that 
the number of interference fringes scales with the phase 
space separation of the classical peaks. This is precisely 




FIG. 1. Wigner function of a coherent state in a localized 
at (a, m a ) = (7r, 10). The function is constant in the other 
angles and a Kronecker delta in the other angular momenta. 



the behavior known from the standard Wigner function 
of superposed Glauber coherent states. 

Molecular alignment dynamics. — The following appli- 
cation aims at demonstrating the viability of the pre- 
sented phase space approach also to describe dynamical 
situations. Specifically, we study the nonadiabatic align- 
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FIG. 2. Wigner function of a superposition of two a-coherent 
states, localized at (a,m a ) — (it, 4) and (n,— 4). The number 
of fringes in the interference contribution is 8, as expected 
from Am a — 8. The blue shaded colors indicate negative 
values (a = 1). 



merit of symmetric top molecules initiated by picosecond 
laser pulses [25]. In these experiments the expectation 
value of cos 2 f3 quantifies the alignment. The dynamics 
after the initial laser kick is governed, in the rigid rotor 
limit, by the field-free Hamiltonian of the symmetric top, 
thH = AS 2 + (C - A) J 2 , where A = H/2h and C = h/2I 3 . 
Using the conventions of [55] , we describe the initial kick 
by a Gaussian pulse at i = with duration 10 and 
strength Q R = 10 (time in units of hjl\). 

Figure [3] shows how the Wigner function of the molec- 
ular orientation evolves in time, taking 1 333) as initial 
state. The half-moon shape in snapshot (a), shortly af- 
ter the kick, is due to the /3-dependence of the interac- 
tion Hamiltonian Hj nt = Qjja cos 2 $ . The molecular 
alignment achieved can be assessed from the form of the 
quasi-probability distribution. One also notices an inter- 
ference structure between the two arms of the half-moon, 
indicating the coherence in the state. The following three 
snapshots (b), (c), and (d) show the Wigner function at 
fractions of the period associated with the classical tra- 
jectory centered on the half-moon. One observes that 
the distribution function by and large follows the classi- 
cal trajectories (gray lines) at these short times, though it 
starts to disperse. In snapshot (e), at a much later time, 
the distribution is spread over the accessible phase space 
and the alignment is lost. Remarkably, at t = n (and in 
multiples of thereof) the alignment revives, marked by a 
small dispersion in /3, see Fig. [3](f ) . One observes from 
the Wigner function that this phenomenon can be asso- 
ciated with an angular momentum-superposition of two 



localized wave packets. 




Ji/4 7C/2 3 71/4 7C/4 7C/2 3 71/4 n 



P P 

FIG. 3. Phase space representation of molecular alignment 
dynamics. The red (blue) shades indicate positive (negative) 
values; gray lines represent classical phase space trajectories. 
Panel (a) depicts the state right after the kick, corresponding 
to maximal alignment. Panels (b), (c), and (d) show how the 
rotational wave packet evolves at the short times i = 0.05, 
0.09 and 0.13. After dispersion, at the later time i — 1.44, 
it is completely delocalized and alignment is lost, see (e). A 
quantum interference effect related to the eigenvalue spacing 
of the symmetric top leads to an alignment revival at t = 
nir, n G N. Panel (f) provides the Wigner function at the 
first revival, while (a) is recovered for i — 2n. (W (SI, m) is 
constant in a, 7 and a Kronecker delta in m a ,m 7 ) 

Conclusions. — In summary, a quantum phase space 
representation has been introduced for the rotation dy- 
namics of rigid bodies. Based on the natural canonical 
variables, i.e. the Euler angles and the associated angular 
momentum projections, it inherits all the relevant prop- 
erties of the standard Wigner function. This includes 
its interpretation as a quasi-probability distribution with 
correct marginals, its agreement with the Weyl corre- 
spondence rule and with the semiclassical limit. Provid- 
ing an appropriate and intuitive framework for discussing 
quantum effects in the rotational motion, we expect it to 
find applications in molecular physics and beyond. 
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